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We show that coexisting domains of coherent and incoherent oscillations can be induced in an 
ensemble of any identical nonlinear dynamical systems using the nonlocal rotational matrix cou¬ 
pling with an asymmetry parameter. Further, chimera is shown to emerge in a wide range of the 
asymmetry parameter in contrast to near ^ values of it employed in the earlier works. We have 
also corroborated our results using the strength of incoherence in the frequency domain {S^j) and 
in the amplitude domain {S) thereby distinguishing the frequency and amplitude chimeras. The 
robust nature of the asymmetry parameter in inducing chimeras in any generic dynamical system 
is established using ensembles of identical Rossler oscillators, Lorenz systems, and Hindmarsh-Rose 
(HR) neurons in their chaotic regimes. 

PACS numbers: 05.45.-a, 05.45.Xt, 89.75.-k 


I. INTRODUCTION 

A fascinating emergent phenomenon in a network of 
nonlocally coupled identical oscillators is characterized 
by the simultaneous existence of synchronous and desyn- 
chronous domains of oscillators [2]. This phenomenon 
was first observed by Kuramoto and Battogtokh in a 
model of complex Ginzburg-Landau equation with expo¬ 
nentially decaying nonlocal coupling m- Later investi¬ 
gations by Abrams and Strogatz revealed similar states in 
an ensemble of phase oscillators with nonlocal coupling, 
who coined the term “chimera” for such states with coex¬ 
isting coherent and incoherent domains [5]. Since then, 
identifying chimera states and their existence criteria has 
become an active area of research both theoretically and 
experimentally [316]. 

In particular, chimera states were observed in networks 
of identical, symmetrically coupled Kuramoto phase os¬ 
cillators, coupled Stuart-Landau oscillators with nonlocal 
interactions in one and two-dimensional arrays and also 
in systems with delay coupling IM]. Recently, it has 
been shown that a symmetry breaking coupling in Stuart- 
Landau oscillators leads to the manifestation of chimera 
death m- Very recently, the existence of chimera 
states have also been demonstrated experimentally in 
populations of coupled chemical oscillators, in electro- 
optical coupled map lattices by liquid-crystal light mod¬ 
ulators [24l [23 and in a mechanical experiment involv¬ 
ing two sub-populations of identical metronomes coupled 
in a hierarchical network [23 • Real world examples ex¬ 
hibiting the states mimicking chimera states include the 
unihemispheric sleep of animals [13, power grids m, 
and so on. Understanding intricacies involved in the dy¬ 
namics of chimera states is also very important from the 
perspective of neuroscience, as it is believed that it can 
be associated with the concept of “bumps” of neuronal 


activity [TS] , 

Initial investigations on the phenomenon of chimera 
have adopted an ensemble of phase-only models in the 
weak coupling limit and demonstrated that chimera 
will arise when the ensemble is distributed spatially us¬ 
ing a nonlocal coupling with a well tuned Sakaguchi 
phase-lag parameter, which introduces a phase asym¬ 
metry Enann]. In particular, it was established that 
chimera states can arise typically when the phase-lag pa¬ 
rameter is near |. Later investigations have extended the 
notion of chimera beyond the phase-only models both 
in the case of nonlocal coupling as well as global cou¬ 
pling [2Qll22] . Interestingly, a couple of recent investi¬ 
gations [131123] have considered dynamical systems that 
can be reduced to phase models and employed a nonlo¬ 
cal rotational matrix coupling with a phase-lag parame¬ 
ter, similar to the Sakaguchi phase-lag parameter. Gor- 
respondingly, the dynamical systems in their periodic 
regimes exhibit transitions, for appropriate value (near 
|) of the phase-lag parameter, observed in the phase- 
only models including chimera states. However, it is not 
clear whether the nonlocal rotational matrix coupling can 
indeed induce chimera in a generic dynamical system that 
may not be reduced to phase models or even in dynami¬ 
cal systems that can be reduced to phase models but in 
the aperiodic (chaotic) regimes. 

Motivated by the above, in this paper we consider the 
collective dynamics of ensembles of identical nonlinear 
oscillators interacting via nonlocal rotational matrix cou¬ 
pling. For vanishing value of the asymmetry (phase-lag) 
parameter in the rotational matrix coupling, the ensem¬ 
ble of dynamical systems is characterized by a direct cou¬ 
pling and it exhibits a global synchronization above a 
threshold value of the coupling strength. In contrast, 
for the ensemble is driven by a cross (conjugate) cou¬ 
pling rendering the ensemble to exhibit asynchronous os- 
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FIG. 1. (Color online) Time averaged frequency of the oscillators (top row) and spatiotemporal plots (bottom row) as a 
function of the oscillator index i — = 500 for an ensemble of Rossler oscillators with nonlocal rotational matrix 

coupling for coupling radius r = 0.3, coupling strength £ = 0.1 and for different values of the asymmetry parameter a. (a) 
& (b) completely asynchronous oscillations for a — 1.57, (c) & (d) coexisting domains of coherent and incoherent oscillations 
illustrating frequency chimera for a = 1.47, (e) & (f) frequencies are synchronized while there exists a drift in the amplitude 
of the oscillations, which are clearly indicated by domains of coherent and incoherent oscillations in the snap shot depicted in 
the inset, for a = 0.97 corresponding to amplitude chimera, (g) & (h) cluster states with frequency entrainment for a = 0.57, 
and (i) & (j) the ensemble of oscillators are completely synchronized for a = 0.1. 


cillations for any value of the coupling strength. For in¬ 
termediate values of the asymmetry parameter, we find 
a wide range of collective dynamical behaviors includ¬ 
ing two different types of chimeras, thereby elucidating 
the emergence of chimeras in a wide range of the asym¬ 
metry parameter in contrast to the earlier studies stick¬ 
ing to near ^ values of the asymmetry parameter. In 
particular, when the asymmetry parameter is increased 
from a zero value, for appropriate coupling strengths, 
we find the occurrence of both frequency and ampli¬ 
tude chimeras in generic dynamical systems that can¬ 
not be reduced to phase models. In addition, a gallery 
of collective dynamical states including clusters, coher¬ 
ence and fully synchronous oscillations arise in these sys¬ 
tems. On the other hand, in the earlier investigations 
frequency chimera alone is shown to exist in the periodic 
regimes with nonlocal rotational matrix coupling in sys¬ 
tems that can be reduced to phase models m- Further, 
the chimera states are characterized by inhomogeneous 
spread of coherent and incoherent domains unlike homo¬ 
geneous spacing of coherent and incoherent domains in 
phase models and in periodic oscillators. We establish 
these facts using three paradigmatic models, namely cou¬ 
pled Rossler oscillators, Lorenz systems, and Hindmarsh- 
Rose (HR) neurons with nonlocal and rotational matrix 
coupling. 

The plan of the article is as follows. An ensemble of 
identical Rossler oscillators will be introduced in Sec. II, 
where the role of the asymmetry parameter in the nonlo¬ 
cal rotational matrix coupling in inducing frequency and 
amplitude chimeras will be discussed. Section III deals 
with the quantitative characterization by introducing the 
strength of incoherence in frequency and amplitude do¬ 
mains thereby distinguishing frequency and amplitude 
chimeras, respectively. Global bifurcation scenario will 



FIG. 2. Snapshots of the variables Xi (left column) and the 
time averaged frequency of the oscillators (right column) of 
the ensemble of Rossler oscillators for the coupling radius 
r = 0.3, the asymmetry parameter a = 1.47 and for different 
values of the coupling strength £. (a) incoherent amplitude for 
£ = 0.05, (b) coherent frequency for £ = 0.05, (c) coexisting 
coherent and incoherent domains of amplitude for e = 0.22, 
(d) coexisting coherent and incoherent domains of frequency 
for £ = 0.22, (e) coherent oscillations with fixed phase differ¬ 
ence with the nearby oscillators for s = 0.35 and (f) coherent 
frequency for £ = 0.35. 


be discussed in Sec. IV and the generic nature of the 
asymmetry parameter in inducing chimeras in a wide 
range of it will be established using ensembles of Lorenz 
oscillators and Hindmarsh-Rose neurons in Sec. V. Fi¬ 
nally a summary of the results and discussions will be 
provided in the concluding section VI. 
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II. FREQUENCY AND AMPLITUDE CHIMERA 
IN AN ENSEMBLE OF ROSSLER OSCILLATORS 

We consider the following ensemble of Rossler oscilla¬ 
tors with nonlocal, and rotational matrix coupling 

^ j=i+P 

Xi=-yi-Zi + — ^ [aii{xj-Xi) + ai2{yj-yi)\ 

j=i — P 

(la) 

^ j=i+P 

yi=Xi + ayi + — ^ [a2i{xj - Xi) + a22{yj - Vi)] 

j=i—P 

(lb) 

Zi = h + {xi - c)zi (Ic) 

where i = 1,2,..., A = 500. The system parameters 
are chosen as a = 0.42, 5 = 2 and c = 4, for which the 
individual Rossler oscillators exhibit chaotic oscillations. 
£ is the coupling strength and P G (1, A/2) is the number 
of nearest neighbors on each side of any oscillator in the 
ring with a coupling radius ^ which provides the 

measure of nonlocal coupling. The coefficients aik^ where 
l,k G 1,2, are the components of the rotational matrix 
coupling represented as 

<^11 <^12 \ _ f cosa sina 
<^21 <^22 J ~ \ —sina cosa 

Here a is the asymmetry parameter. The matrix B facil¬ 
itates both the direct coupling as well as the cross cou¬ 
pling between the variables x and y. For a = one 
has the cross (conjugate) coupling where the y variable 
is coupled to x and xtoy and the ensemble is in complete 
asynchrony (see Figs, [^a) and (b)) in the entire range 
of 5. On the otherhand, the coupling is direct for a = 0 
with the X variable coupled to x and ^ to for which 
the ensemble of coupled systems, Eq. 0, is in complete 
synchrony for 5 > 0.1 (see Figs, [^i) and (j)), while for 
intermediate values of a, different other interesting dy¬ 
namical states occur (see Figs.[^c) and (h)) as described 
below. Thus the asymmetry parameter a plays a crucial 
role in determining the collective dynamics of the ensem¬ 
ble of coupled systems with nonlocal rotational matrix 
coupling. For intermediate values of the asymmetry pa¬ 
rameter a G [0, |], the coupling in Eq. Q includes both 
direct and cross couplings with the value of the weights 
aik determining the proportion of them thereby dictat¬ 
ing the collective dynamics of the ensemble. For suit¬ 
able intermediate values of the asymmetry parameter, 
the ensemble exhibits a mixture of the above extremi¬ 
ties thereby leading to the emergence of coexisting coher¬ 
ent and incoherent domains by a spontaneous splitting of 
the ensemble. In particular, coherence and incoherence 
can be found both in the amplitude and frequency, and 
also incoherence in the amplitude alone but with coher¬ 
ent frequency, thereby distinguishing frequency chimera 
and amplitude chimera, respectively. Considering chaot¬ 
ically evolving coupled oscillators, the frequency chimera 
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FIG. 3. Same as Fig. [^but for the value of the asymmetry 
parameter o = 0.95. In this case, the time averaged frequen¬ 
cies of the oscillators are in coherence, for all values of the 
coupling strength we have examined, as shown in the insets 
of the figures, (a) asynchronous oscillations of the ensemble 
of Rossler oscillators for £ = 0.03, (b) coexisting domains of 
synchronous and asynchronous oscillations for e = 0.1, (c) co¬ 
herent oscillations for e = 0.2, and (d) complete synchronous 
oscillations for e — 0.3. 


is characterized by both temporal and spatial chaos in 
the incoherent interval whereas only spatial chaos is ob¬ 
served in the incoherent interval of the amplitude chimera 
while the corresponding temporal dynamics being peri¬ 
odic in most cases m- Beside the above two types of 
chimeras, cluster states, coherent states and completely 
synchronized states are also observed during the dynam¬ 
ical transitions as discussed below. 

We have fixed the coupling radius as r = 0.3 to inves¬ 
tigate the dynamical transitions by gradually decreasing 
the asymmetry parameter from a = ^ = 1.57 for the 
value of coupling strength 5 = 0.1. The ensemble of 
identical oscillators exhibits a complete asynchrony for 
Of = I as discussed above for any value of the coupling 
strength. The cross (conjugate) coupling induces a drift 
in the oscillations of identical oscillators resulting in the 
distribution of their frequencies in a large range as shown 
in Figs.j^a) and (b). The values of uji for each oscillator 
is calculated as uJi = 27rKi/At^ i = 1, 2,3,..., A, where 
Ki is the number of maxima of the time series Xi{t) of 
the oscillator during the time interval At. Decreasing 
the asymmetry parameter to a = 1.47, a small fraction of 
the direct coupling competes with the conjugate coupling 
in reinforcing certain oscillators to oscillate in synchrony 
resulting in a group of synchronized oscillators while re¬ 
taining the rest in asynchrony (see Figs, [^c) and (d)). 
This scenario of spontaneous splitting of the ensemble 
into coexisting coherent and incoherent domains is noth¬ 
ing but the chimera state. It is evident from Figs, [^c) 
and (d) that both the frequency and the amplitude of 
all the oscillators are in synchrony in the coherent do¬ 
main, while in the other domain they remain desynchro¬ 
nized thereby characterizing the underlying state as a fre¬ 
quency chimera. It is also clear from these figures that 
there is an inhomogenity in the spread of the coherent 
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and incoherent domains unlike the homogeneous spread 
of both the domains in phase only models and in the 
case of periodic oscillators. Upon decreasing a further 
to (a = 0.97, the proportion of the direct and the conju¬ 
gate couplings are close to each other and in this region 
the ensemble is entrained to a single frequency resulting 
in complete coherence (see Fig. [^e)). Nevertheless, the 
amplitudes of all the oscillators are not completely locked 
resulting in partially coherent and incoherent domains in 
amplitude as depicted in the inset of Fig. [^e), which 
is also clearly visualized in the corresponding spatiotem- 
poral plot in Fig. [^f). This state with the coexistence 
of oscillators with coherent and incoherent amplitudes 
while their frequencies are completely entrained is clas¬ 
sified as an amplitude chimera. Synchronized clusters 
emerge for further decrease in the asymmetry parameter 
as illustrated in Figs, [^g) and (h) for a = 0.57, where 
the fraction of the direct coupling is higher than that of 
the conjugate coupling. When the direct coupling dom¬ 
inates the conjugate coupling for further lower values of 
the asymmetry parameter the ensemble achieves com¬ 
plete synchronization, which is shown in Figs, [^i) and 
(j) for a = 0.1. 


Now we will trace the dynamical transitions by vary¬ 
ing the coupling strength for a couple of fixed values of 
the asymmetry parameter a. Snapshots and average fre¬ 
quencies of the ensemble of oscillators are shown in the 
left and right column of Fig.[^ respectively, as a function 
of the oscillator index for a = 1.47 for different values of 
the coupling strength. The amplitude of the oscillators 
are completely uncorrelated for e = 0.05 as can be seen 
in Fig. [2](^, even while their frequencies are entrained 
(see Fig. ^b)). Upon increasing the coupling strength 
further, frequency chimera comes into existence as char¬ 
acterized by coexisting asynchronous and synchronous 
domains in both the amplitude (Fig. [^c)) and the fre¬ 
quency (Fig. I^d)) of the oscillators. Finally, all the os¬ 
cillators evolve in coherence as depicted in Fig. [^e) with 
synchronized frequency (see Fig. [^f)). Next, we will fix 
the asymmetry parameter at a = 0.95, where we observe 
the amplitude chimera. Snapshots along with the fre¬ 
quency profile of all the oscillators are depicted in Figs.j^ 
The frequencies of the oscillator ensemble are always en¬ 
trained for the chosen value of the asymmetry parameter 
as can be seen from the insets of Figs. The oscillators 
evolve independently for 5 = 0.03, which is evident from 
the random values acquired by the oscillator ensemble 
depicted in Fig.J^a). Amplitude chimera emerges for 
5 = 0.1 (see Fig.l^b)), where the oscillators are grouped 
into synchronous and asynchronous domains depending 
upon the amplitude correlation. Eventually, all the os¬ 
cillators evolve in coherence with fixed phase difference 
with the nearby oscillators for further higher values of 
the coupling strength, the snapshot of which is shown 
in Fig. [^c) for 5 = 0.2. Finally all the oscillators are 
in complete synchronization for 5 > 0.23 as depicted in 
Fig. [^d) for 5 = 0.3. 


III. QUANTITATIVE CHARACTERIZATION 
OF CHIMERA STATES IN THE ENSEMBLE OF 
ROSSLER OSCILLATORS 


Recently, we have introduced the strength of incoher¬ 
ence S [21] , as a measure to distinguish various collective 
dynamical states, which is defined as 


5 ' = 1 - 


s 

2^m=l 

M 


Sm = 0{S - (Ji{m)), (3) 


where ©(•) is the Heaviside step function, and 6 is a pre¬ 
defined threshold. Here, we take as a certain percent¬ 
age value of difference between the upper/lower bounds, 
^i,i,maxlmin^ of fhe allowed values of xi/s. M = 20 is 
the number of bins of equal size n = N/M = 25, where 
N is the total number of oscillators. The local standard 
deviation is introduced as 


^ mn 

= (, - E [zi,j-< rn = 1,2,.. .M, 

\ j=n{m—l)-\-l 


(4) 


where Zi^i = Xi^i — I = l,2...(i, d is the dimen¬ 

sion of the individual unit in the ensemble, i = 1,2...V, 
< zi^m >= and (...)( denotes av- 

erage over time. When is less than 5^ Sm = 1, 

otherwise it is ‘O’. In the incoherent domain, the lo¬ 
cal standard deviation cr/(m) has some finite value Vm, 
which is always greater than the predefined threshold 6 
and hence Sm = 0,Vm, thereby resulting in unit value 
for the strength of incoherence S in this domain. On 
the other hand, in the case of coherent domain the stan¬ 
dard deviation ai{m) is always zero in the coherent do¬ 
main and hence Sm = l,Vm, thereby resulting in a zero 
value for S for this domain. Since the chimera states 
are characterized by coexisting coherent and incoherent 
domains, the strength of incoherence S will have interme¬ 
diate values between zero and one, 0 < S' < 1. However, 
the strength of incoherence is incapable of distinguish¬ 
ing different types of chimera. In order to facilitate this, 
we estimate the strength of incoherence in the frequency 
domain S^;, as different from S where we use the same 
prescriptions as above with the choice zi^i = uji^i — 

Now, can be used to clearly distinguish amplitude and 
frequency chimera. Since the frequencies of all the oscil¬ 
lators in the ensemble is entrained for amplitude chimera, 
the strength of incoherence in the frequency domain for 
the amplitude chimera is always 5'^^ = 0, whereas it varies 
between 0 < 5'^^ < 1 for frequency chimera. 

S and Soj are depicted in Figs, [^a) and (b) as a func¬ 
tion of the coupling strength 5 for two different values 
of the asymmetry parameter a = 1.47 and 0.95, re¬ 
spectively, characterizing the dynamical transitions dis¬ 
cussed in Figs. and for the specific coupling radius 
r = 0.3. In the range 5 G (0, 0.09) the unit value of S 
indicates that the oscillators evolve independently while 
at the same time = 0 corroborates that the frequency 
of the oscillators are synchronized. Fluctuation in the 
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FIG. 4. (Color online) (a) and (b) Strength of incoherence in amplitude S and in frequency Suj as a function of the coupling 
strength characterizing the dynamical transitions discussed in Figs. [^and|^ respectively, for r = 0.3. Time averaged frequency 
of the oscillators are shown in the insets. In (c) and (d), S and Su are plotted as a function of the coupling radius r for the 
fixed coupling strengths e = 0.2 and 0.1, respectively. The value of the asymmetry parameter has been fixed as a = 1.47 in (a) 
and (c), while it is chosen as a = 0.95 in (b) and (d). Degree of disorder in the oscillator frequencies is increased in the range 
of frequency chimera as evident from the insets in (a) and (c). 




FIG. 5. (Color online) Plots of the spatial correlation func¬ 
tion G(0) — G{i) (a) for amplitude chimera for r = 0.3 and 
a = 0.95, showing a power-law dependence, and (b) for fre¬ 
quency chimera for r = 0.3 and a = 1.47, showing no power- 
law dependence. 


values of S and S^j between null value and unity con¬ 
firms the existence of frequency chimera in the range of 
5 G (0.09,0.33). Beyond 5 = 0.33 both S and S^j attain 
null value corresponding to the existence of complete co¬ 


herence among the oscillators. It is to be noted that the 
frequency profile of the oscillators is distributed over a 
wide range in the frequency chimera state thereby show¬ 
ing an increase in the degree of disorder among the oscil¬ 
lator frequencies (see the inset of Fig. [^a)). As discussed 
earlier in Fig. [^for a = 0.95, the frequency of the oscil¬ 
lator ensemble is always entrained as indicated by the 
null value of Suj in Fig. Qb), whereas up to 5 = 0.03, 
S acquires unit value attributing to asynchronous ampli¬ 
tude variation. Amplitude chimera exists in the range 
5 G (0.03,0.16), where S fluctuates between zero and 
unity while = 0 in this range. All the oscillators 
evolve in coherence for 5 > 0.16, which is confirmed by 
the null values of S and Suj in Fig. [^b). 

Next, the strengths of incoherence S and are shown 
in Figs, [^c) and (d) as a function of the coupling radius 
r for two different values of the asymmetry parameter 
a = 1.47 and 0.95, respectively. The intermediate values 
of S and S^j between zero and one in Fig. I^c) for 5 = 0.2 
characterize the existence of frequency chimera in a wide 
range of the coupling radius r G (0.07,0.49). It is to be 
noted that for r = 0.5, the nonlocal coupling turns out 
to be the global coupling. The value of S^j in Fig. |^d) 
for e = 0.1 remains zero while S fluctuates between zero 
and unity in the range of r G (0.09, 0.49), elucidating the 
existence of amplitude chimera as a function of the cou¬ 
pling radius. Thus, it is clearly evident from both these 
figures that chimera states can indeed appear in a rather 
wide range of the coupling radius r. Further, the fre¬ 
quency profile of the oscillators is distributed over a wide 
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chaos, whereas it does not have a power-law distribu¬ 
tion in the case of frequency chimera (see Fig. [^b)). 
We have also checked that spatial correlation function 
indeed exhibits a power-law distribution in the entire pa¬ 
rameter space of the amplitude chimera corresponding 
to Fig. 4(b). Nevertheless, we have confirmed that both 
the amplitude and frequency chimeras are characterized 
by a discontinuous peak at the origin corroborating the 
existence of spatial incoherence as demonstrated by Ku- 
ramoto and Nakao. 


0 0.25 0.5 

e IV. GLOBAL BIFURCATION DIAGRAM FOR 

THE ROSSLER OSCILLATORS 



FIG. 6. (Color online) Two parameter phase diagram for the 
dynamical states of the Rossler oscillators 0 as a function of 
the coupling strength e and the asymmetry parameter a. See 
Table I and text for details. Note that the power-law scaling 
turbulence exists only in the amplitude chimera region (AC) 
as explained in the text. 


TABLE I. (Online color) Characterization of different collec¬ 
tive states with the corresponding color coding in Figs. 6,8 
and 10: 


Dynamical state 



So 


Remarks 

DSYC 

1 

0 

1 

1 


FC 

c 

c 

c 

d 

0 < c, d < 1 

AC 

c 

0 

c 

d 


CL 

c 

0 

0 

d 


CO 

0 

0 

0 

d 


CS 

0 

0 

0 

0 



range in the frequency chimera state thereby showing an 
increase in the degree of disorder among the oscillator 
frequencies as shown in the inset of Fig. [^c). 

In order to check how the amplitude and frequency 
chimeras characterized by the strengths of incoher¬ 
ence as discussed above relate to the power-law scal¬ 
ing turbulence in the nonlocally coupled chaotic sys¬ 
tems studied by Kuramoto and Nakao in Ref. [29], we 
have estimated the spatial correlation function G{i) = 

= 0> 1> •••f- “ = Xj, where N 

is the total number of oscillators in the ensemble. It 
has been demonstrated in [29] that the power-law spatial 
correlations appear generically in self-oscillatory media 
with nonlocal coupling exhibiting spatiotemporal chaos. 
In particular, for small but finite i, the correlation func¬ 
tion G{i) behaves like G( 0 ) — G{i) « 7o 0 < <^ < 2 , 

where 70 and 71 are constants, exhibiting a power-law 
dependence on the distance and a discontinuous peak 
(70 7 ^ 0 ) at the origin when the coupling constant de¬ 
creases below a critical value. We have found that the 
spatial correlation function G{0) — G{i) exhibits a power- 
law dependence on the distance i as shown in Fig. HJa) 
for the case of amplitude chimera characterizing spatial 


The global scenario including the discussed dynam¬ 
ical transitions as a function of the coupling strength 
in the range 5 G (0,0.5) and the asymmetry parameter 
a G (0,1^) is depicted in Fig. The parameter regions 
for which the oscillators are in the desynchronized state, 
cluster state, coherent state, and complete synchroniza¬ 
tion are marked as ^DSYG\ ‘GL’ ‘GO’ and ‘G5'’, re¬ 
spectively, while that for the amplitude and frequency 
chimeras they are indicated by ‘AG’ and ‘FG’, respec¬ 
tively. The corresponding color codeing and the nature 
of various strengths of incoherence are explained in Table 
I. 

In Figj^ the parameter space is demarcated using the 
strength of incoherence S to distinguish the incoherent, 
coherent and amplitude chimera regimes. is used to 
distinguish the frequency chimera from the amplitude 
chimera. Method of removal of discontinuity [21] is used 
to distinguish cluster states from the other states. In 
this case, we introduce another measure Sq which is the 
strength of incoherence after removal of discontinuities, 
while the previous measure S corresponds to the strength 
of incoherence before the removal of discontinuity. The 
values of S and So indicates the existence of cluster state. 
The strength of incoherence Sx obtained by using the val¬ 
ues of the dynamical variables xi^i directly instead of the 
difference variables zi^i = xi^i — distinguishes be¬ 

tween complete synchronization (GS) with Sx = 0, F = 0 
and coherent state (GO) with Sx = S = 0, 0 < d < 1. 
We have used the same set of initial conditions, which is 
uniformly distributed between 1 and —1, for all combi¬ 
nations of the parameters in the entire parameter space. 

For a = I, the ensemble of Rossler oscillators evolves 
asynchronously for the entire range of the coupling 
strength. In the range of the asymmetry parameter a = 
(1.55,0.97) transition from desynchronized state to fre¬ 
quency chimera and finally to coherent state is observed 
as the coupling strength is varied in the range 5 G (0, 0.5). 
Decreasing a further, thereby increasing the proportion 
of direct coupling in the range a = (0.97,0.6), the en¬ 
semble of oscillators transit to amplitude chimera from 
asynchronous state as the coupling strength is increased 
from zero. Further increase in the value of 5 leads to 
complete coherence among the oscillators and finally to 
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FIG. 7. (Color online) (a) Space time plot and (b) time aver¬ 
aged frequency as a function of the oscillator index for s = 3.5, 
and (c) S and Su as a function of the coupling strength elu¬ 
cidating the frequency chimera in an ensemble of identical 
Lorenz oscillators with nonlocal rotational matrix coupling 
for r = 0.35 and a = 1.46. 
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FIG. 8. (Color online) Two parameter phase diagram for the 
dynamical states of the Lorenz oscillators 0 as a function of 
the coupling strength e and the asymmetry parameter a. See 
Table I and text for details. 


complete synchronous oscillations. For 0 < a < 0.6, we 
observe similar dynamical transitions except for the fact 
that the amplitude chimera is surrounded by the cluster 
states in both the directions of the coupling strength as 
seen in Fig. 


V. CHIMERA IN ENSEMBLES OF LORENZ 
OSCILLATORS AND HINDMARSH-ROSE 
NEURONS 


In this section, we will establish our claim that the 
rotational matrix coupling with the asymmetry param¬ 
eter a is capable of inducing chimera in any dynamical 
system for a wide range of values of a by demonstrating 
the existence of chimera in ensembles of identical Lorenz 
oscillators [36] 


^ j=i+P 

Xi = < 7 {yi - Xi) + — ^ [au{xj-Xi) + ai2{yj-Vi)], 
j=i-P 
^ j=i+P 

yi=Xi{p-Zi)-yi +— ^ [a 2 i{xj - Xi) + a 22 {yj - Vi)], 

j=i-P 

Zi=Xiyi-bzi, ( 5 ) 

and Hindmarsh-Rose (HR) neurons [37] 

Xi=yi- axi +bx‘f - Zi + I ( 6 ) 

^ j=i+P 

+ ^ H - Xi) + ai2{yj - yi)], 

j=i-P 

^ j=i+P 

yi=c-dxl-yi +— ^ [a2i{xj-Xi)-\-a22{yj-yi)\, 

j=i-P 

Zi =ro[s{xi - xo) - Zi], ( 7 ) 

where i = 1,2,..., A = 500. In Eq. ([^ a = 10, 
p = 28, 6 = 8/3 are the system parameters of the in¬ 
dividual Lorenz oscillators, while in Eq. 0a = c=l, 
d = 5, ro = 0.005, s = 4, xq = —1.6, 5 = 3, / = 3 
are the parameters of the HR neurons. Eor the chosen 
parameter values, the individual Lorenz oscillators and 
HR neurons exhibit chaotic behavior. The spatiotem- 
poral plot and the average frequencies of the individual 
oscillators are shown in Eig. [^a) and (b), respectively, 
for a = 1.46, r = 0.35 and 5 = 3.5 elucidating the ex¬ 
istence of frequency chimera in the ensemble of Lorenz 
oscillators. It is to be noted that the frequency chimera 
is indeed a multichimera with four coherent and inco¬ 
herent domains as can be clearly visualized in both the 
spatiotemporal and the average frequency plots. To be 
more clear, the oscillators simply populate the two dif¬ 
ferent lobes of the attractor resulting in four coherent 
domains (see Eig. IZla)), as can be readily seen in the 
time averaged frequency uui in Eig. IZlb), separated by 
the incoherent borders due to non-local coupling. Note 
that chimeras and multichimeras can also be further dis¬ 
tinguished quantitatively, by introducing the measure of 
discontinuity as shown in [21] , though this is not pursued 
here. The strength of the incoherence S for the ampli¬ 
tude and in the frequency domain for the ensemble of 
Lorenz oscillators are depicted in Eig.j^c) as a function of 
the coupling strength e G (0,5). The null(unit) value of 
Suj{S) confirms that the frequencies(amplitudes) are syn¬ 
chronized (desynchronized) in the range 5 G (0,1.9). The 
fluctuation in the values of the strength of incoherence S^j 
and S corroborates the existence of frequency chimera in 
the range of the coupling strength 5 G (1.9,4.6). Eurther, 
the inset of Eig. [^c) indicates an increase in the ran¬ 
domness in the frequencies of the identical Lorenz oscil¬ 
lators for the values of 5 where we observe the frequency 
chimera. Beyond 5 = 4.6, the zero values for both and 
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FIG. 10. (Color online) Two parameter phase diagram for 
the dynamical states of the HR neurons 0 as a function of 
the coupling strength e and the asymmetry parameter a. See 
Table I and text for details. 


FIG. 9. (Color online) (a) Space time plot and (b) time aver¬ 
aged frequency as a function of the oscillator index for £ = 0.3, 
and (c) S and Soj as a function of the coupling strength elu¬ 
cidating the frequency chimera in an ensemble of identical 
Hindmarsh-Rose (HR) neurons with nonlocal rotational ma¬ 
trix coupling for r = 0.35 and a = 1.46. 

S indicate the emergence of complete coherence among 
the oscillators of the ensemble. 

A two parameter bifurcation diagram as a function 
of the coupling strength 5 G (0,5) and the asymme¬ 
try parameter a G (0, |) of the ensemble of identical 
Lorenz oscillators is shown in Fig. We have used 
the same symbols as in Fig. to indicate the differ¬ 
ent emergent behaviors. As in the case of the ensem¬ 
ble of Rossler oscillators, the ensemble of Lorenz oscil¬ 
lators evolve asynchronously in the entire range of the 
coupling strength e G (0,5) ioi a = ^. In the range 
of the asymmetry parameter a G (1.57,1.05), there is a 
transition from asynchronization (DSYC) to frequency 
chimera (FC) and then to coherent states {CO) followed 
by complete synchronization (CS) as a function of the 
coupling strength. Above a = 1.05, the dynamical tran¬ 
sition is in the sequence of desynchronization (DSYC)^ 
amplitude chimera (AC), cluster states {CL), coherent 
states {CO) and complete synchronization {CS) upon 
increasing the coupling strength from null value. 

The existence of frequency chimera (multichimera), 
that is coexistence of coherence and incoherence domains 
in the ensemble of HR neurons for a = 1.46, r = 0.35 and 
5 = 0.3 is shown in the spatiotemporal plot (see Fig. ia)) 
and the average frequency plot (see Fig. gb)). The 
dynamical transition from asynchronous state to com¬ 
plete coherence via frequency chimera in the ensemble 
of HR neurons is evident from the values of S^j and S 
in Fig. I^c). The coherent and incoherent oscillations 
in frequency and amplitude, respectively, are revealed 
from the values of 5'^ = 0 and S' = 1 in the range 
5 G (0,0.09). Intermediate values of S^j and S, that is 
between null value and unity, in the range of the cou¬ 
pling strength 5 G (0.09,0.51) confirm the existence of 


frequency chimera. As in the ensemble of Rossler and 
Lorenz oscillators, the frequency chimera is accompanied 
by an increase in the randomness in the frequencies of the 
ensemble of identical HR neurons. The ensemble of HR 
neurons are in complete harmony for 5 > 0.51 as is clear 
from the values of S^j and S from Fig.j^c). A two param¬ 
eter phase diagram as a function of the coupling strength 
5 G (0,0.6) and the asymmetry parameter a G (0, |) of 
the ensemble of HR neurons is shown in Fig. [OT We 
have followed the same abbreviations as in Figs.|^and[^ 
to indicate the different dynamical regimes here also. All 
the HR neurons evolve independently for a = | in the 
entire range of the e we have analyzed. There is a tran¬ 
sition from asynchronous state to synchronous state via 
frequency chimera followed by coherent states as a func¬ 
tion of the coupling strength for a G (1.57,0.86). For 
further larger values of the asymmetry parameter, the 
frequency chimera is replaced by the amplitude chimera 
and cluster states in the above dynamical transition as 
seen in Fig. 

VI. CONCLUSION 

In conclusion, we have shown that the nonlocal rota¬ 
tional matrix coupling can indeed induce chimera states 
in an ensemble of a generic nonlinear dynamical systems 
exhibiting chaotic oscillations that cannot be reduced to 
phase models, whereas dynamical systems in the periodic 
regimes facilitating phase reduction have been employed 
so far along with such coupling. In particular, using an 
ensemble of identical Rossler systems, exhibiting chaotic 
oscillations for appropriate parameter values, we have 
demonstrated that frequency and amplitude chimeras are 
induced by the nonlocal rotational matrix coupling in a 
rather large range of the asymmetry parameter a as a 
function of the coupling strength in contrast to the ear¬ 
lier studies limiting to near | values of it in the phase- 
only models and in the dynamical systems that can be 
phase reduced. Further the chimera states are charac- 
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terized by inhomogeneous spread of coherent and inco¬ 
herent domains unlike homogeneous spread of these do¬ 
mains in phase models and in periodic oscillators. In 
addition, a range of collective behaviors including cluster 
states, coherent states and complete synchronization are 
observed during the dynamical transitions as a function 
of the system parameters. We have also corroborated 
the emergence of frequency and amplitude chimeras us¬ 
ing the strength of incoherence in the frequency domain 
(Suj) and in the amplitude domain (S'). The generic na¬ 
ture of the nonlocal rotational matrix coupling in induc¬ 
ing chimera in any nonlinear dynamical systems in a wide 
range of the asymmetry parameter is also confirmed by 
illustrating the emergence of chimeras in an ensemble of 
identical Lorenz systems and Hindmarsh-Rose neurons 
in their chaotic regimes. Further, increase in the ran¬ 
domness of the frequencies of the ensemble of identical 
oscillators may be a signature of the frequency chimera 
as observed in all the three paradigmatic models we have 
analyzed. Thus it is clear that the asymmetry parameter 
in the nonlocal rotational matrix coupling, which induces 


the asymmetry in the ensemble of identical nonlinear sys¬ 
tems, plays a crucial role in determining the desired dy¬ 
namical behavior. Hence for appropriate choice of the 
asymmetry parameter, it is possible to obtain the desired 
collective states in dynamical systems including neuronal 
systems [34], hydro dynamical systems such as laminar 
and turbulent regions in Couette flow studies [32] l33] . 
etc., with appropriate nonlocal coupling. 
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